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ABSTRACT 

We present four-color images of CRL 2688 obtained in 2009 using the Wide- 
Field Camera 3 on HST. The F606W image is compared with archival images 
in very similar filters to monitor the proper motions of nebular structure. We 
find that the bright N-S lobes have expanded uniformly by 2.5% and that the 
ensemble of rings has translated radially by 0'.'07 in 6.65 y. The rings were ejected 
every lOOy for ~4 millennia until the lobes formed 250y ago. Starlight scattered 
from the edges of the dark E-W dust lane is coincident with extant H2 images 
and leading tips of eight pairs of CO outflows. We interpret this as evidence that 
fingers lie within geometrically opposite cones of opening angles ~ 30° like those 
in CRL618. By combining our results of the rings with ^^CO absorption from 
the extended AGB wind we ascertain that the rings were ejected at ~18 km 
with very little variation and that the distance to CRL 2688, Vexp/dexp, is 300- 
350 pc. Our 2009 imaging program included filters that span 0.6 to 1.6//m. We 
constructed a two-dimensional dust scattering model of stellar radiation through 
CRL 2688 that successfully reproduces the details of the nebular geometry, its 
integrated spectral energy distribution, and nearly all of its color variations. The 
model implies that the optical opacity of the lobes > 1, the dust particle density 
in the rings decreases as radius"^ and that the mass and momentum of the AGB 
winds and their rings have increased over time. 

Subject headings: stars: AGB and post-AGB, stars: winds, outflows, (ISM:) planetary 
nebulae: individual (CRL 2688), (ISM:) reflection nebulae 
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Introduction 



CRL 268^is the ideal pre-planetary nebula ("pPN") for studying the historical ejection 
of mass from many AGB and post-AGB stars. Optically, its bright ~ 15" nebular core 
contains a polar (N-S) pair of bright expanding lobes and a dark equatorial (E-W) dust 
lane, both of which are presumed to have formed at about the same time as the transition 
from the AGB to post-AGB states of late stellar evolution. These features are characteristic 
of pPNe in the DUPLEX subclass (functionally, those that contain deep equatorial dust 



extinction) (Ueta et al. 2000 Sahai et al. 2011 ^ In addition an ensemble of concentric 



quasi-circular arcs (hereafter "rings") of total radial extent > 40" surrounds the lobes and 
the dust lane. The arcs are spaced semi regularly with characteristic separations ~ 1". The 
high brightness of the lobes and the rings as well as their proximity to make CRL 2688 a 
particular productive site to study their characteristics. 

CRL 2688 has a wealth of complementary observational information. For example, the 
visible and near-IR illumination of CRL 2688 consists of direct and scattered light from a 
carbon- and nitrogen-rich but completely obscured supergiant of apparent spectral type F5 



la with Teff ^ 6500 K (Klochkova et al. 2000 Cohen & Kuhi 1977 Crampton et al. 1975). 



Lying adjacent to the nebular center of symmetry are smaller structures seen in the radio 



continuum (Jura et al. 2000), the mid infrared (Le Mignant et al. 2007; Goto et al. 2002 



hereafter "G02"), and a cavity of diameter 2" seen in ^^CO (2-1) (Cox et al. 2000, hereafter 
"COO"). Like the star, these regions are fully obscured within the dust lane at near IR 
and visible wavelengths. (We largely ignore them in this paper since they do not appear 



1 aka Egg Nebula, RAFGL2688, V* V1610 Cyg; RA2000 = 21 02 18.8, ^2000 = +36 41 41^ 
l,b = 080.2-06.5 

^ see also the compilation of images of pPNe at 
pPNe/ 



www . astro . Washington . edu/balick/ 
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to affect the structure of HST images.) Further afield, the bases of both of the lobes are 
especially prominent in the thermal infrared (G02), suggesting that starlight scatters into 
the lobes through dense and warmed dust that might affect its initial colors. Finally, H2 
emission outlines the perimeter of the dust lane and the tips and edges of the N-S lobes (Le 



Mignant et al.||2007| [Kastner et al.||2001[ |Sahai et al.||1998a[ hereafter "S98a"). The shocked 
H2 traces the interaction zones between the (presumably) expanding lobes and dust lane, 
on the one hand, and the ambient halo containing the rings, on the other. 

Lying outside of the lobes and the dust lane, the well illuminated, delineated, and 



resolved rings of CRL 2688 represent the early history of AGB mass ejection. Hrivnak et al 



(2001) found similar but less conspicuous rings in IRAS 16594-04656 and IRAS20028+3910. 



Other rings have been seen in IRAS 17441-2411 by Su et al. (1998), IRAS17150-3224 by 



Kwok et al.| ( |1998D , IRAS19114+0002 by Su 200^ IRC+10216 by [Mauron fc Huggins 



(2000). Of these, the rings of CRL 2688 are the best illuminated and defined. However, 



aside from their density modulations little is known of the ages, expansion patterns, and 
ejection histories of any of these ensembles of rings. 

The broad focus in this paper is the evolution patterns of the lobes and the rings. 
Thus we briefly review observations of their kinematics here. So far the only single-dish 



CO studies can detect the flux from the entire ensemble of the rings (e.g., Bujarrabal et al. 



2001, hereafter BOl, & Hrivnak & Bieging 2005). Obviously such studies cannot resolve 



the expansion pattern. Spatially resolved kinematic studies of CRL 2688 consist primarily 
of ^^CO observations. Cox et al. (2000, hereafter "COO") found several pairs of fingers 



confined within the lobes and the dust lane. (COO elected to associate each of the fingers 
with coUimated "jets".) However, they could not detect the low-brightness halo with the 



www. astro .Washington. edu/users/balick/APN/APN_tmeeting_tdocuments/ 



Tuesday/ksu . ppt 
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rings in their study. A few years later Fong et al. (2006) observed ^^CO in CRL2688 using 



closer antenna spacings and better sensitivity to large structures. They marginally detected 
the inner portions of the halo containing the rings, but could not fit an expansion pattern 
to this zone. 

One of our primary science objectives is to probe the structure of mass loss at the 
times when the lobes and rings first formed through measurements of the proper motions of 
these features. This science objective is plausible because, as we assert next, the structures 
in CRL 2688 seen today accurately mimic their geometry at the time of their formation. 
First, we note that pPNe do not show signs of extended nebular emission lines from species 
such as (though a few show low-ionization emission lines from localized shocks). This is 
a crucial observation since the onset of nebular ionization drives strong D-type ionization 
fronts into the neutral gas that are preceded by an expanding supersonic "piston" that can 
quickly sweep up, displace, and shock heat the upstream gas. 

Additionally, unlike mature PNe, there are no traces of rapidly expanding central 
cavities that are formed by fast winds (~ 10'^ km s~^) that create a hot, sparse central 
"bubble" (10^ — lO^K). The supersonic growth of the bubble displaces and compresses the 



gas that once occupied its interior (e.g., Kwok et al. 1978, Schonberner et al. 2005, and 
subsequent papers). Also, there is no evidence of any other forces (e.g., strong embedded 
magnetic fields) or substantial shock-induced radiative energy losses that impede or alter 
the structures that have been observed in pPNe. Therefore it comes as no surprise that 



the dense, bright cores of pPNe are found to expand uniformly (Alcolea et al. 2001 



Castro-Carrizo et al. 2002 Sanchez Contreras et al. 2004 & Alcolea et al. 2007). In other 



words, the initial geometries and ejection speeds of the lobes and rings of CRL 2688 are 
preserved in the refiection nebulae observed now. 



There are very few previous studies of proper motions in pPNe and none of their rings. 



(Borkowski & Harrington 2001) were the first to study proper motions in pPNe. They 
used a pair of images separated by ~ 3 y to monitor the apparent motions of knots near 
the symmetry axis of Hen 3-1475. CRL2688 has no counterparts of these knots. Of more 



direct relevance, Ueta et aL (2006, hereafter "U06") used a pair of HST NICMOS images of 
CRL2688 at 2fim separated by 5.5 y (unspecified camera and spatial resolution) to study 
proper motions of the N-S lobes. U06 derived an inclination i of 7.7° from the plane of the 
sky, a distance D of 420 pc, a stellar luminosity L of 10^'^ Lq, a total shell mass M of 1.2 
M0, and an expansion age Tkin of 350}]^ The north lobe is tipped forward. At this distance 
1" corresponds to 6.3 xlO^° km = 0.0020 pc. 

CRL2688 is close enough (see below), many of its structures are sufficiently bright and 
sharp, and its foreground extinction is small enough (Ay = 0.88 mag) that if a feature 
has an in-sky motion of 10 - 20 km s~^ then changes in the structure of its rings can be 
monitored at high spatial resolution with telescopes such as the Hubble Space Telescope 
("HST") over 5 — 10 y. In this study we used high-resolution images from 2002 and 2009 
at 0.6yum from the HST archives. Compared to the NICMOS images used by U06, our 
data have a smaller diffraction core (FWHM ^ 0'.'067 after dithering vs > 0'.'2, depending 
on NICMOS camera), lower readout noise and dark current (i.e. better sensitivity to the 
faint rings), and a slightly longer time base (6.65 y). Indeed, we successfully measured the 
aggregate proper motions of sharp features in the rings as well as the lobes. The images are 
described in §2, analyzed in §3, and combined with complementary ^^CO data in §5. 

A second focus of this paper is an analysis of the scattered light of CRL 2688 as 
observed in four colors, 0.6, 0.8, 1.1, and 1.6/im in 2009. These images are used to constrain 
a new and detailed model of radiative transfer by dust in order (1) to infer properties of 



^ Note that other inchnations and distances 16° and 1 kpc, respectively, have been sug- 
gested by lYusef-Zadeh et al.l (|1984D ) . 
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the dust and its distribution and (2) to fit the spectral energy distribution ("SED") of 
CRL2688. See §4. Tliere liave been two prior attempts of radiative transfer modeling of 



CRL2688 (Skinner et al. 1997; Lopez & Perrin 2000). Skinner et al. used a dusty outflow 



with hollow polar lobes, but their geometry was not constrained by high resolution imaging 



that emerged shortly thereafter (S98a, Sahai et al. 1998b, hereafter "S98b", and Weintraub 



et al. 2000). Lopez & Perrin proposed a model based on a dust density law smoothly 



decreasing with latitude from the midplane. However, deep images of CRL 2688 do not 
support this geometry. 

Our major findings are summarized in §6 and a brief discussion of viable ring formation 
mechanisms is presented. 



2. Data Calibration 

We obtained dithered images of CRL 2688 using WFC3 (Wide Field Camera 3 
Instrument Handbook, Version 2.1, Dressel et al. 2010) through four broad-band filters, 
F606W, F814W, FllOW, and F160W (Table 1.) Detailed characteristics of the camera 
and its filters are described on the WFC3 web page^ We obtained dithered short- and 
long-exposure images through each filter. The F606W and F814W exposures were limited 
to a single HST orbit in August 2009. FllOW and F160W exposures were obtained during 
another orbit in October 2009. All exposures were made with a 512x512 subarray in order 
to expedite data transmission to the ground. The visible images were initially drizzled onto 
a 0'.'0396 grid, the native pixel size of the UVIS camera on WFC3. A proper motion of 10 
km s~^ per decade corresponds to the width of a full F606W-F814W pixel at the assumed 
distance of CRL 2688 of 420pc. The infrared images were drizzled to a grid of 0'.'128. A few 

J 

iwww. stsci . edu/hst/wf c3 
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saturated pixels in some images of the lobes were replaced with unsaturated values obtained 
from short exposures. 

Color overlay: The color distribution of scattered light was measured from flux- 
calibrated WFC3 images in each filter after using bright field stars to align each image. 
The flux calibration is based solely on calibration data supplied with the images from the 
standard calibration pipeline. The basic results are shown in Fig[T| In this figure the color 
saturation and contrast were altered in order to exaggerate the brightness of the faint outer 
rings and to emphasize the overall color variations. 

Proper-motion studies: Only images in the F606W filter were used for measuring 
proper motions since the archival data in this filter are the highest in spatial resolution 
and signal-to-noise ratios. Table 2 shows a log of the F606W images including exposures 
dates and duration and the original pixel sizes and fields of view. Detecting changes 
in structure based on images originally made using different cameras at various epochs 
requires thoughtful image alignment and gridding, the best possible corrections for camera 
sensitivity variations and geometric image distortions, careful image registration, and 
well-matched point-spread functions. We describe our procedures next. 

In every case we used the latest corrections supplied by STScI staff in May 2010. 
Note that geometric distortion in the central 512^ subfields of the UVIS and IR cameras 
are very small. Distortions in the ACS/HRC camera have been carefully measured since 
2002. Corrections for distortion were applied to the pipeline data from the archives using 
Multidrizzle (described in the January 2009 version of the handbook]^ Based on HST 
documentation we believe that the systematic errors in the corrections are generally much 
less than 0.1 pixel (we have no means of confirming this result). In short, large-scale spatial 

^ www. stsci . eciu/hst/HST_overview/documents/multidrizzle/multidrizzle .pdf 
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distortions that might materially affect our results should be insignificant. 

Image alignment is tricky since no suitable stationary point-like sources are found 
in the field of the CRL 2688 images. That is, each of the field stars in the images have 
idiosyncratic motions relative to CRL 2688. As noted earlier, the nucleus of CRL 2688 
is not visible within any of our filter passbands. Thus the images from each epoch and 
filter were initially aligned from the absolute positions of the images that are stored in 
the image headers. (These coordinates are based on guide stars, each with its own proper 
motion.) Then image pairs were aligned empirically using the centers of symmetries of 
large-scale features, such as the lobes and arcs. We are confident that the uncertainty of the 
image alignment, ~ C'DOS, is a small fraction of the PSF of the cameras (FWHM = 0'.'067 
after dithering). The nebular symmetry center agrees closely with the polarization center 



(Weintraub et al. 2000) and the ^^CO kinematic center (COO) 



We searched for proper motions of the structures within CRL 2688 in all image pairs. 
All of our pair-by-pair results and our assessment of their uncertainties are compiled in 
Table 3. Note that CRL 2688 was observed using the low-resolution (O'.'l/pixel) WF3 chip 
in WFPC2 in which the PSF varies noticeably with chip location. Unfortunately, the 1995 
image is not sufficiently deep for reliable results of changes in the positions of the faint 
outer rings. Thus the result for the 1995-2009 image pairs in Table 3 is based largely on 
the innermost rings. What's more the camera focus appeared to be poor since we had to 
convolve the 2009 images to a PSF of 0'.'2 to match the PSFs of field stars in the 1995 image. 
Using the 1998 and 2009 image pairs produced rather somewhat irregular inconsistent 
patterns of nebular ring growth. In short, our assessment of the systematic errors led us to 
concerns about the quality of the results obtained using any of the early WFPC2 images for 
studies of large-scale proper motions that our error estimates may not capture. 

The results derived from ACS/HRC (c. 2002.94) and WFC3 images (c. 2009.59) are 
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the most credible and the only ones discussed hereafter (6.65y time base). Those images 
were first drizzled onto a 0'.'025 grid (the native pixel size of the ACS/HRC image) and the 
intensities were normalized to the fluxes of three bright field stars. As a check, we confirmed 
that this image renormalization resulted in an image ratio that is very close to unity in the 
smooth, bright portions of the lobes, as expected. Therefore any residuals from multi-epoch 
image subtractions are not artifacts of errors in camera flux calibrations. 

Finally, the 2002 and 2009 images were converted from counts per second to units of 
energy flux, divided (after suitable magnifications to correct for nebular expansion), and 
examined for changes in surface brightness. The conversion to energy flux is based on the 
value of PHOTFLAM stored in the image header. PHOTFLAM accounts for the system 
and filter throughput and the widths of the filters. Values of PHOTFLAM are based on 
methodical calibrations of the system by STScI staff and are very stable for the F606W 
filters used in our study. 

3. Analysis & Results 

3.1. Image structures and color variations 

The structure and color variations of CRL 2688 were presented in Fig [1} The most 
prominent new features are these: 

• the southern lobe and the pairs of "searchlights" that continue beyond its tips are 
noticeably redder than the northern lobe and its corresponding pair of searchlights. 

• the opening angle of the searchlights widens and the gaps in their centers close 
noticeably at longer wavelengths. We surmise that the gap is formed by a shadowing 
obstacle whose edges become more transparent with increasing wavelength. 
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• the rings are generally bluer in color than the lobes. One exception is where the rings 
cross the southern lobe that is tipped backwards. 

• the dark dust lane and especially the finger tips at its extremities are outlined in 
continuum light in our l.l/im image. This outline is far less prominent in our 0.8/im 
and 1.6/im images. 

• the north (south) lobe faded by 10-12% (2-4%) between 2002 and 2009. The E-W 
dust faded with an irregular pattern. The northern (southern) quadrant of the lobes 
faded by 15-20% (5-15%). See Fig 2b and 2c for the pattern of brightness changes 
between 2002 and 2009. 



N-S lobes: The lobes are fiame-shaped emitters that both absorb and scatter starlight 
radiated into their base. As such, they are probably the main source of illumination for the 
surrounding ensemble of rings at lower latitudes. 

The north lobe is marbled with deep extinction filaments (most with H2 emission 
counterparts). Its relatively blue colors are partially attributable to the shorter interior line 



of sight to its forward-tipped lobe. Also, as pointed out by Cohen & Kuhi (1980), the north 
(south) lobe is seen in forward- (back-) scattered light. It will appear bluer if the scattering 
phase function is not isotropic. 

The fainter southern lobe is generally smoother in appearance than its northern 
counterpart and devoid of marbled dust. It divides cleanly into two mirror-symmetric 
features separated by a gap along its symmetry axis that extends into the searchlights. We 
surmise that the south lobe is hollow, or that the interior gap along its symmetry axis is 
an illumination shadow cast from the E-W dust lane, or both. The continuous dark gap 
extending from the base of the southern lobe into the searchlights support the shadow 
hypothesis without eliminating the idea of a hollow lobe core. It is difficult to tell whether 
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the north lobe is similar in this regard. 

Searchlights: The bright searchlights are best explained if the lobes are semi-transparent 
to starlight when viewed from very high latitude. The red gaps in the searchlights are 
shadows, as proposed by S98b and others. Each searchlight has a dark core along the 



nebular symmetry axis. As proposed by S98b and Le Mignant et al. (2007), H2 knots at the 
tips of each lobe may be occulting disks that cast shadows into each searchlight. In practice, 
the exact location of the symmetric pair of dark occulting disks cannot be located in any of 
our images. They could be anywhere on the symmetry axis, perhaps even somewhere inside 
the dark E-W dust lane. 

E-W dust lane: We noted that the edges of the dust lane are outlined nicely in our 
1.1/im— and, to a lesser extent, in the 0.8 and 1.6 — /im images. The scattered light follows 
essentially the same brightness pattern as the shock-excited H2. As noted above, the 
fingertips are especially bright. Given the faintness of the rings at low latitude, it seems 
unlikely that the fingertips are illuminated from the exterior. 

Close to the geometric center of nebular symmetry the present infrared images show 
two compact and very reddened features that have no visible counterparts. Both of these 
lie near the outer edge of the central CO cavity of 2" in diameter. The brighter of these IR 
features was identified as "knot A" by G02 in their studies at 3.8 and 4.7/im. The same 
knot appears in Fig [T] and was seen at 1.65 and 2.15/im by S98a. The fainter and redder of 
the two near-nuclear infrared features aligns with a feature at the base of the southern lobe 
reported by G02. These features are also very proximate to the presumed locations of the 
two "polar caps" in the scattering model to be described in §4. 
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3.2. Growth Patterns of the Lobes 

Methodology: Monitoring changes in the locations of sharp structures must be tailored 
to the nebular geometry; e.g., changes in a circular feature involve finding a change in 
its radius. Our efforts to fit the lobe expansion rely on the sharpness of the edges of the 
lobes and the filamentary dust features in their interior. The earlier image of a registered 
pair was magnified a factor (1 + M), where M > is a free parameter. The ratio of the 
magnified earlier image to the newer image of the pair produces a "ratio image" in which 
radial growth is seen as large-scale dark-light features wherever sharp gradients in intensity 
are located^] 

The signal of structure change is clear when no magnification is applied (M = 0). M 
was incremented until the systematic signals of change were minimized (i.e., "fiattened"). 
The goodness of fit is measured by eye since standard statistical measures of variations are 
dominated by readout noise, hot pixels, and the strong residuals left by translated field 
stars. We found that for the eye, a small range of values of M fiattens the residuals equally 
well. This range determines the uncertainty of the fit. 

Results: The residual images are shown in Fig. [2] for three values of M applied to 
the ACS/HRC image of 2002.94 prior to forming the ratio image: (a) no magnification 
(M = 0), (b) (M = 0.009) and (c) (M = 0.025). In the case of the lobes the residuals in the 
ratio image are very nicely fiattened for M = 0.025 ± 15%. In other words, A6r oc Or, where 
A6r is the radial offset of a lobe feature at a radial offset form the symmetry center at Or = 

^ As a rule, we found that the image ratios are a better measure of alignment than are the 
differences, at least for the outermost rings. This is because the surface brightness falls off 
rapidly with radius,. So the pattern of differences at the outer edges of the images become 
very faint. However, the rings can be traced to the noisy outer edges of a ratio image. 
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0. (N.B.: We were unable to tell whether the same value of M applies to the expansion of 
the edges of the E-W dust lane.) 

This analysis yields an expansion age for the lobes of 6.65y/0.025, or 266y ±15% in 
good agreement with the estimate of (1) 350y found by U06 based on their comparison 



of lower-resolution HST/NICMOS images spanning 5.5y, (2) 200y found by Jura & Kroto 



(1990) for the time since the star left the AGB (assuming a distance of 1 kpc), and (3) 
the CO kinematic ages of 52 to 500 y for various pairs of CO knots embedded in the lobes 
(COO after adjusting to D = 420 pc). Earlier age estimates require a prior knowledge of the 
inclination angle i and a distance D and, hence, are less accurate 



3.3. Growth Patterns of the Rings 

Methodology: We investigated two models for fitting the growth pattern of the rings: 
(1) a brief epoch of ring ejection followed by a homologous "uniform expansion" in which 
Or oc 9], and (2) "ongoing ring ejections" with constant motions in which Or oc 9^ for each 
ring. The former suggested itself during image alignment, as noted in §2. Uniform expansion 
is the most common Doppler expansion pattern in many other pPNe and PNe (§1). On 
the other hand, the ongoing ring ejection is a popular but thus-far untested paradigm in 
which an AGB star is somehow regularly stimulated to eject a bubble of enhanced density 



® Ages derived from proper motions are not affected by uncertainties in i or D. However, 
age estimates based on measured Doppler shifts are very sensitive to any errors in D and 
sinz, especially when i is small. U06 estimated a lobe inclination oi i = 7.7° but could not 
rule out larger values. All prior estimates of i and D vary from those of U06 up to a factor of 



two (Yusef-Zadeh et al. 1984). We note that i ~ 15° is nicely consistent with our scattering 



model in 54. 
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or a new compression wave every 10^ — lO^y (§1) at the same speed (so that the rings do 
not overtake one another and lose their separate identities). As noted in §1, this idea is 
motivated by the ensembles of concentric rings seen in many other pPNe and PNe. 

In order to find the expansion pattern of the rings, the observed images were 
transformed into polar coordinates and their alignments carefully checked before forming 
the image ratio. As seen in Fig. [3} circular rings and arcs become horizontal lines. On the 
other hand, if the rings each have the same ejection speed Vj. = DO,, then the entire system 
of rings will translate by a fixed amount A6r in the radial direction between observing 
epochs. Thus the residual will be best fiattened by applying an optimum choice A6r for the 
radial translation of the ring ensemble of the earlier image prior to forming the image ratio. 

Results: Both models of uniform expansion fit the WFC3 data well. We find that the 
ratio image is optimally fiattened everywhere if we adopt M = 0.009 ± 0.002. This value 
of M corresponds a the kinematic age of the ensemble of rings is Atfcj„ = 6.65y/M = 740y 
with 20% uncertainty. A very good fit is also found for the constant ring ejection-speed 
model. A translation A9r = O'.'OT ± O'.'Ol nulls most of the ensemble as seen in Fig. [sj 
However, strong residual structure is quite conspicuous on both sides of the dark E-W dust 
lane within the white box shown in this figure. (Ignoring this region presumes that the 
structure within the white box is part of a different kinematic system such as a lobe or the 
adjacent dark dust lane.) ^^CO observations discussed in §5 will break this degeneracy in 
the quahty of the models. 

As noted in §1, a study of the growth pattern of the rings is a key result of this paper. 
However, the present study is limited to arcs within ^ 20" of the nucleus. Therefore 
an analysis of multi-epoch images spanning the full ~ 80" field of the rings seen in the 
ACS/WFC F606W image of 2002 would be a very useful next step for a better description 
of the ensemble expansion pattern and the ejection histories of the rings of CRL 2688. 
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4. Model for Scattered Light 

4.1. Overview and Previous Results 

The most conspicuous features of the reflection nebula CRL 2688 seen in HST images 
are its highly articulated and symmetric shape. The WFC3 images trace only scattered 
starlight, some of it multiply scattered. Clearly this provides a highly biased view of the 
actual density distribution. Therefore a better assessment of the spatial distribution of dust 
requires a realistic analysis in which starlight computationally transferred through a series 
of dusty features of different geometries. In this section we describe a successful model in 
which each structural component is generally characterized by the geometric features in 
our images and the dust is described by particles of various sizes, size distributions, and 
scattering and absorption properties. Details of the model can be found in Appendix A. 



4.2. Model Methodology, Strategy, and Assumptions 



Constructing a model that simultaneously reproduces the complex structure, color 
variations, and SED of CRL 2688 is challenging. We have employed the radiation transfer 



code LELUYA (Vinkovic 2003, see also www.leluya.org). LELUYA is a powerful 
yet flexible tool which nicely solves the integral equation of the formal solution of 
radiative transfer including dust scattering, absorption and thermal emission in two 
spatial dimensions. The transfer equations are solved on a highly unstructured triangular 
self-adaptive grid that simultaneously traces both the density gradients and the optical 
depth gradients over many orders of magnitude in spatial and optical depth space. It allows 
the user to assign grains of various sizes and particle distributions to each feature. However, 
LELUYA does not treat polarization and non-isotropic scattering. 



We use spherical amorphous carbon dust grains (Hanner 1988). The optical constants 
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for these grains are derived by presuming isotropic scattering and Mie theory. For simphcity 
we assume that the chemical properties of the dust particles are uniform; i.e., only their 
size distribution varies from feature to feature. 

Constrained and unconstrained model variables: Fig. |4] provides a conceptual sketch 
of the structural components used for the computations. We tried to specify as many 
modeling constraints from the observed four-color images (§2) as possible in order to 
reduce degeneracy of the results of the model. However, the core of the nebula where the 
illumination originates is hidden from direct view even beyond 4/im (G02) so its structure 
was initially assumed and then optimized to fit the SED and the patterns of scattered light. 
We next describe how the many model variables were fixed or fitted. 

The dust distribution is shown schematically on the left side of Fig |4] and described 
in detail in the Appendix. A total opening angle of 30° was adopted for the lobes and the 
searchlights, 10° for the shadow zone along its symmetry axis, and 40° for the polar caps. 
The density profile adopted for the caps, the ellipsoidal equatorial disk, the hollow lobes, 
and the envelope scales as r°, r~°'^, r"*"^, and r~^, respectively. 

The geometry and density of the polar caps and the radiating area of the star were 
varied in order to match the SED and color distributions of the structures. We adopted 
two models for the illumination source with different temperatures, 4000 K and 7000 K. 
Scattering and absorption of starlight was then followed through the dust features. The 



particle-size-dependent opacities are taken from Manner (1988) 



Other features are those identified in S98a and the circumstellar cavity noted in COO. 
The polar caps enclose this cavity close to the star where dust particles initially form. The 
stellar illumination is launched through annular openings in the caps on its way to the lobes 
and the searchlights. A distance 420 pc was assumed (§1). We treated z as a free parameter 
to be determined by assuming that the two lobes are of equal brightness and given the dust 
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distribution in the model envelope, tilting the lobes to match their observed brightness 
ratio. 

It is worth noting that we are unable to derive even an approximate dust mass for 
several reasons. One is that the grain sizes (specifically, the scattering area per unit dust 
mass) are not tightly constrained by the model. Also, we assume a uniform grain radius or 
size distribution for each structural component. This is obviously unrealistic. Moreover, 



we do not try to incorporate scattering by large dust particles into the model. (Jura et al 



2000) estimated that the large grains could be at least as large as 5 mm. Such large grains 



may have a major impact on the derived dust mass. 



4.3. Model Results 

The success of the model in fitting the data is shown in Figs. |4] (right side) and [s] 
Generally speaking, the overall correspondence of the geometry and color variations of the 
model and our HST images is excellent. Of course, much of this is the direct result of the 
flexibility in designing the density and dust particle size distribution of the model. The 
resulting color variations are in good agreement except (arguably) in the faintest outer 
rings and searchlights where the quality of the images is poorest. The modeled arcs show 
the basic morphological properties of imaged arcs, even though we use ideal spherical shells 



in our model. We confirm the estimates by Sahai et al. (1998b) that arcs represent a local 



density increase by a factor of 3. A key feature of the model is that the SED is reproduced 
very successfully across the spectrum from the visible to the sub-mm. 

We note that the long run time of the model made a detailed parameter study of many 
parameters impossible. Thus it is difficult to assess the uniqueness and the uncertainties on 
the selection of all of the parameters in the model. 
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N-S lobes and E-W dust lane: We find that the optical opacity r of the lobes ~ 2 
at 0.55/im and that the opacity of the disk is ~55 in its midplane. The model image is 
shown at an inclination angle i = 15° . This value of i best compensates for the unequal 
observed surface brightness ratio of the north and south lobes (assuming that both lobes 
are identical and that scattering is isotropic). The result agrees with the value of i = 16° 



derived from the scattering model for CRL 2688 by Yusef-Zadeh et al. (1984). Even so, the 
fit to the image morphology and SED is adequate for 5° < 2 < 15°. A list of other derived 
model parameters is provided in able |4] of Appendix A. 

Central star: We obtained a bolometric flux of 1.0 ± 0.2 x 10"^ W m~^. This result is 
based on a correction for asymmetric diffuse radiation. The correction arises from the lack 
of spherical symmetry which is typically ignored in the literature when the luminosity is 
calculated. That is, more energy is visible from a vantage point along the symmetry axis 
(large inclination angles) than near the midplane where the disk is blocking the central star 
(small inclination angles). Our model takes this into consideration. A detailed discussion of 



this effect is given by Vinkovic et al. (2004) in their §A1. 



We estimate the light originates from a star with a luminosity of 5500 ± 1100 L© and a 



temperature T^, ~ 7000 K, in close agreement with Klochkova et al. (2000) assuming D = 
420 pc and Ay = 0.88. The corresponding stellar radius is 0.23 ± 0.02 A.U. (51 ± 5Rq), or 
0.56 msec ±10%. 

Circumnuclear core: The radiative transfer calculation yields the radius of the dust-free 
cavity of 46.4 stellar radii, deep within the CO cavity. This value critically depends on the 
dust condensation temperature, which is in our case Tjnax=1000K in the midplane. Since 
this value is not a measurement, and because models with other temperatures are viable, 
we adopted this value to reduce model degeneracy. Newly formed dust radiates in the near 
IR spectrum which is completely obscured by the dust lane. Thus photometric verification 
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of the size, temperature, and radiated flux of the is not amendable to empirical verification. 

To first order, the computed dust temperature distribution is dominated by a simple 
gradient outside of the polar caps near the star. The temperatures in all directions 
drop smoothly from 1000 K at log(r/i?i„) ~ to 300K at log(r/i?i„) ~ 1 and lOOK at 
log(r/_Rj„) ~ 2. Here is the dust condensation radius (the radius of the inner grid). The 
dust temperature continues to fall smoothly to the edge of the entire grid at Rout = 600R 



With this cavity size and assumed limit on the condensation temperature we can 



rule out the maximum dust temperature of 330K as suggested by Lopez & Perrin (2000). 
Asymmetric diffuse radiation and differences in grain sizes create variations in the highest 
dust temperature on the spherical surface of dust-free cavity. The large 0.4/im grains in 
disk and lobes experience temperatures from lOOOK at the midplane to ~1100K on the 
polar axis. In contrast, the small grains (0.01 - 0.05/xm) in the polar caps reach 1400K 
since smaller grains are relatively inefficient emitters of infrared radiation. 

Polar Caps: Recall that knot A seen in the infrared (G02 and Fig [T]) lies within 1" 
of the expected location of the north polar cap. It corresponds with the bright zone at 
the upper edge of the polar cap. Similarly, the fainter infrared knot lies near the expected 
location of the heavily obscured polar cap at the base of the southern lobe. 

Rings: Although the radiative transfer model was designed to reproduce the geometry 
of the rings, searchlights, and the shadow zone on the nebular symmetry axis, it does not 
fully capture the color variations outside of the nebular core. The HST images indicate that 
the searchlights are relatively blue along their full extent whereas the model shows the blue 
colors are found only where the nearby blue lobes dominate the illumination of the rings. 
Some of this can be attributed to the two-dimensional nature of the model which cannot 
accurately treat multiple scatterings. 
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Two other areas of model- image disagreement are (1) the bright parts of the searchhghts 
and their relatively red edges and (2) the opening angle of the shadow within the searchlight 
that decreases as the wavelength increases. These problems would probably be resolved 
if we model the annular aperture that defines the outer edges of both the lobes and the 
shadow as having "soft" edges with opacity gradients. 

Another area of disagreement is that the disk used in the model is darker and thicker 
than its counterpart seen in HST images. This is readily accommodated by a change in 
the assumed geometry and density of the disk. In reality, the dark disk seems to be highly 
structured (COO, G02 and §3), so full agreement with a highly idealized model agreement is 
not anticipated. 

SED: The agreement of the predicted and model SED is one of the major successes 
of this model. The best fit to the SED is shown in Fig. [5) The fit starts to fall below 
the data at mm wavelengths where thermal radiation from large dust particles dominates. 



Nonetheless, the overall quality of the fit is far superior to those of previous models (Skinner 



et al. 1997; Lopez & Perrin 2000). 



Dust particles: The dust particle size and size distribution were treated as free 
parameters and adjusted to provide a satisfactory match to the observed image data and 
the SED. The data do not tightly constrain these particle sizes. We adopted a uniform 
grain size of OAfim for the disk and lobes, 0.2/im in the envelope including the rings, and a 
particle size distribution that scales as (grain radius) in the polar caps. Available imaging 
and photometry at A > 2^m provides a range of constrains on dust properties, optical 
depths and guidance on the limits of the interpretation of our modeling results. The data 



for the SED were derived from observed fiuxes compiled from Ney et al. (1975) and Skinner 



et al. (1997) and then dereddened using a fixed foreground extinction of Ay = 0.88 (G06). 



Lopez & Perrin (2000) proposed a different radiation-transfer model to explain the 
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features of CRL 2688 based on a dust density law smoothly decreasing with latitude above 
(and below) the midplane. Their modeled images showed brightness spikes in place of 
searchlight beams, which is an interesting result considering the simple configuration of their 
dust geometry. However, these spikes and the overall nebular geometry are simply not in 
agreement with the searchhghts seen in HST images. Even so, one of their firm conclusions 
is that the highest dust temperature is about 330K, far below the dust condensation 
temperature and inconsistent with the nebular SED. 

5. Comparison to CO Distribution and Kinematics 

In this section we combine our results with radio CO studies in order to develop 
additional information on the history of mass ejection in CRL 2688. As it turns out, many 
of the CO, the H2, and the scattered-light features of the lobes and dust lane nebulae agree 
in location nicely, at least in outline. Thus we can expect that outside of the circumnuclear 
CO cavity (see below) the volumes that generate the CO hues and the dust that scatters 
starlight are one and the same. It should be noted that the present HST images cannot 
be used to probe the structure deep within the lobes and the dust lane since the optical 
opacity exceeds unity (§4). 

5.1. Structure and Motions in the Lobes and Dust Lane 

We first consider the lobes and the dark dust lane. In §1 we noted that COO investigated 
both the structure and the kinematics of ^^CO within the lobes and the dust lane of 
CRL 2688 using a radio interferometer with 1" spatial resolution. They found seven pairs of 
finger-shaped outflows bounded by the edges of the lobes and the dust lane. As noted in 
§1, BOl decomposed the Doppler-broadened, nebula-integrated ^^CO and ^^CO line profiles 
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into an unsaturated "fast wind" of width ~ ±40 km s~^ arising in the inner core where the 
finger-shaped CO outflows are found plus a saturated "slow wind" of width ~ ±15 km s~^ 
from the surrounding AGB wind where the rings are located. They and COO also found a 
deep and narrow ^^CO absorption feature whose velocity is blueshifted by —18 km s^^ from 
the systemic velocity of CRK 2688. 

We summarize the high-resolution observations of the lobes and the dark dust lane by 
COO. Seven pairs of fingers are revealed, each pair of which consists of outflows with nearly 
equal and opposite Doppler shifts, as noted in §1. As shown in their Fig. 2, finger pairs 
D1-D2 and F1-F2 lie along the leading edges on the N-S lobe. E1-E2 lie at the intersection 
of the N-S lobes and the E-W dust lane. A1-A2, B1-B2, C1-C2, and G1-G2 are embedded 
within the dust lane. COO suspected that each finger may increase in velocity along its 
length. Bright knots of H2 and scattered starlight at l.l/xm are also coincident at the 
extremities of all of these fingers. 

Interestingly all eight of the finger tips associated with the dust lane also have 
scattered-light counterparts at 0.8 and l.lfim. Some of the edges of the cone are illuminated 
by scattered light from the lobes and the rings in our HST images. However, the outermost 
edges of the of the equatorial dust lane are especially prominent. They are not likely to be 
illumined from the exterior since the rings are particularly faint at low latitudes. So it's 
possible that starlight propagates through hollow CO fingers where it scatters from dusty 
clumps at the fingertips. Speculatively, the fingers (possibly outflowing jets or bullets) lie 
within a cone whose opening angle is ~ 30° - strikingly similar to the morphology of the 



ensemble of finger pairs in CRL618 analyzed by Sanchez Contreras et al. (2002). In other 
words, the equatorial dust ellipsoid adopted for the two-dimensional scattering model in §4 
may instead be a pair of cones with a bundle of three pairs of CO fingers inside. 

COO derived the kinematic ages of each finger pair from the distance traversed from 



- 24 - 



the nucleus and the Doppler velocity: T^^n — '^finger'^^ finger — (D ^rSini)-i-(v/jop/cosi), 
where d finger is the distance traversed by each fingertip, y finger is the corresponding speed, 
D is the distance, and ydop is the measured Doppler shift of relevant tip. They adopted D 
= 1 kpc and i = 16°. If each finger pair has a different value of i, as we suspect, and if D 
= 420 pc, then their estimates of T^j^ need to be revised substantially. Perhaps the best 
CO measure of the kinematic age of the lobes come from finger Gl at the tip of the N lobe 
which, after a correction for D — 420pc, is roughly 500y. In contrast, our estimate of 250y 
of the expansion age for the outer boundary of the lobes is based on proper motions for 
which D and i are not needed. 

COO did not attempt to derive the kinematic age for the dark dust lane. However, the 
separations 6r and Doppler velocities ^oop are very similar for finger pairs A, B, C, and G in 
the dust lane, on the one hand, and D, E, and F in the lobes, on the other. It is noteworthy 
that the fingers, the lobes, and the dust lane may have similar ages in the range of 250 - 
500 y. 

COO also noted the presence of a small shell-hke structure in close proximity to the 
nucleus whose radius is 1". The CO surface brightness drops precipitously within this shell. 
We feel that the CO core could alternately be described as a 1" cavity at the center of the 
core of the CO brightness distribution. The radius of the hole corresponds to the volume 
within which the dust temperature drops from ~ lOOOK at its condensation radius of 26 
msec to about lOOK at falOO times further out (§4). 

5.2. The Kinematics of the Rings 

Without short interferometer spacings COO could not detect the low-brightness, ~ 40" 
halo outside of the core of CRL2688 in which all of the rings are situated. However, 
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Fong et al. (2006) used the BIMA array in conjunction with the NRAO 12-m singe-dish 
telescope to map CRL 2688 in in the ^^CO J = 1 — with ~ 3" resolution. They detected 
CO from extending from an inner radius of 9" out to the noise level at almost 20". Little 
velocity structure in this expanding region was reported. This favors a const ant- velocity 
outflow since uniform expansion would produce a large spread of Doppler shifts. But the 
interferometer's poor sensitivity to very extended structure leaves a large margin of doubt 
in this conclusion. 

BOl, COO, and Fong et al. noted the presence of the ^^CO absorption feature that 
covers the bright core and the more extended halo. As noted first byBOl the absorption 
arises in the extended slow- velocity component. Presumably this absorption region shares 
the same volume as the rings. If so, then the rings lie in a foreground region of low velocity 
dispersion, ±2 km s~^. 

The simplest interpretation of this result is that the entire ensemble rings results 
from repeated ejections, each of them characterized by outflow velocities blueshifted by 
—18 km s~^. The small CO absorption line width implies little variation in outflow speed 
over 4000y. This result is inconsistent with the possibility of brief ejection and uniform 
expansion {9r oc 9l) in §3.3. 



5.3. The Distance to CRL 2688 

If the ensemble of rings of CRL 2688 expand isotropically, then their proper motion 
speed, 9r = O'.'O? ± O'.'Ol in 6.65 y, is the same as the Doppler shift of the ^^CO absorption 
feature from the nebular systemic velocity, = 18 ± 2 km s^^. From this we derive D = 
340 ± 60 pc. The uncertainty is dominated by the measurement of 6r- Alternately, if we 
adopt the width of the slow velocity component of BOl, 15 km s~^, then D = 280 pc. 
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6. Conclusions 

Adopting a distance D = 420 pc and extinction Ay = 0.88 mag, then optical/IR 
continuum, CO, and H2 images can be interpreted consistently as follows (in order of 
increasing distance from the nucleus): 

• The integrated SED and the optical-near IR color distribution is well-fitted with a 
realistic and generally tightly constrained dust density distribution illuminated by a 
cool (~7000 K) and luminous (5500 Lq) star. 

• The central shell (possibly a hole) in CO is associated with very warm dust at its 
condensation temperature needed to fit the nebular structure and the SED. 

• The lobes and dust lane are at least partially filled with finger-shaped outflows. 

• Inclination- and distance-independent expansion ages from the boundary of the lobes 
implies an expansion age fti250y. Based solely on CO kinematics, the fingers within 
the dark dust lane may have a similar outflow age. 

• The polar inclination angle lies between 5 and 15°. We favor larger values in order to 
account for the lower brightness and redder colors of the south receding lobe. 

• The system of rings has a kinematic signature that argues for repeated ejections at a 
constant, regulated velocity, 15-18 km s~^, along the line of sight. 

• Ring ejection continued for several thousand years preceding the formation of lobes. 

• A fit of the scattered starlight in the rings implies a density that decreases as radius^. 

• As a result, the mass and momentum per ring have increased roughly linearly - not 
unexpected as the luminosity and radius of the star increase during the ascent of the 
AGB. 



-27- 



• If the rings expand isotropically then we derive a distance of 340pc ± 60 pc based 
on the measured proper motions of the rings and blueshift of the foreground ^^CO 
absorption feature. 

• At this distance the expansion time between rings - that is their ejection period - is 
lOOy. 

• The entire nebula has faded irregularly by an average of ~ 10% over the course of 
6.65 y. The northern parts of the nebula faded considerably more than the southern. 
The changes are more likely the result of changes in the immediate environment of 
the star than the star itself. 



6.1. Evaluating Ring Formation Scenarios 



Ramstedt et al. (2008) has provided a thoughtful review of the mechanisms that drive 



AGB winds. In addition, they successfully fitted models of constant AGB mass loss Magb 
to the CO profiles of many types of evolved outflow sources. In contrast we find that 
Magb increases (~ linearly) in time. The same result was derived from a careful analysis 



of CO line profile measurements of other pPNe by Hrivnak & Bieging (2005). Thus the 



fundamental issue of the evolution of Magb remains generally unresolved. 

The process for forming the rings in the AGB winds is a matter of considerable 



discussion. As noted by Balick et al. (2001) and others, the rings are too closely spaced 
in time (10^ — 10'^)y, to be the result of thermal pulses since the typical interpulse time is 



10^ - 10 V (cf. Karakas 2010 Karakas fc Lattanzio 2007 for M < 4 Mq and Z = 0.02). The 



periods of surface pulsations are far too short (1-10 y) to produce the arcs. 

The rings must form at the stellar surface or at the base of the flow since the opacity 
of the wind to driving radiation rapidly decreases with radius. Their properties were 
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predicted with the appropriate period using a two-fluid hydro model in which recurring 
large-scale dynamic instabilities alter the grain properties in the winds of a single AGB 



star (Simis et al. 2001). This concept was extended to two dimensions and applied to C 



stars by Woitke (2006). Woitke argues that the extra degree of freedom generates local flow 



instabilities that produce multiple mushroom-shaped ejection events, the leading edges of 



which produce observable arcs. Soker (2000) proposed that the rings form as the result of 



periodic cycles of AGB magnetic activity that modulate the effective surface temperature 
of the star and, hence, the dust formation rate. 

Alternately, the rings have been attributed to the periodic effects of a close companion 
star in an eccentric orbit with a small periastron and the appropriate period to explain 



the separations of the rings. Harpaz et al. (1997) introduced a class of models in which a 



companion star or planet passes inside the grain formation region (very close to the star) 
so as to interrupt the grain formation and, hence, to modulate the density of the wind. 



He (2007) and Raga et al. (2011) proposed that if the close companion has an orbital 



eccentricity e < 0.8 then a spiral wind outflow is extruded. When the observer lies near the 
equatorial plane, as in CRL 2688, the projected outflow pattern morphs to a series of round 
arcs on opposite sides of the that plane that alternate in radial displacement by even or 
odd multiples of the ring half-separation (that is, phased by 180°). Insofar as CRL 2688 
is concerned, the arcs are not sufficiently regular to test this hypothesis rigorously (Fig. Isl). 



Finally, Zijlstra & Bedding (2002) pointed out that the time scales of period changes 



in long-period "meandering Miras" (~ 15% of all Miras) and the inter-ejection times of 
rings are similar. They propose that agreement in these times scales implies "Mass loss 
[sic] variations induced by period changes may be one of the causes of the rings " and they 



^ J. Cernicharo (private communication) has proposed a similar model for the rings of 
lRC+10216. 
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cite R Hya as an illustration. (No physical mechanism that connects the phenomena is 
suggested in that paper.) While the exposed surface of the central star of CRL/,2688 is 
presently too hot to be a long-period Mira, it may have been as recently as a few centuries 
ago when the lobes were ejected. 
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Fig. 1.— Color overlays of the F606W-F814W (left) and F110W-F160W images (right) 
shown to the same scale and pixel size (O'.'OS) with a highly non-linear contrast function 
and exaggerated color saturation. The native resolutions of the images are 0'.'040 (visible) 
and (y.'128 (IR). The FOV is 400 pixels or 32". The color saturation has been enhanced to 
highlight the color changes. 
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Fig. 2.— Ratios of flux-calibrated HST images of CRL2688 in the F606W filter in 2002.94 
(ACS/HRC) and 2009.59 (WFC3). The intensity display scale ranges from 0.5 (black) to 2.0 
(white) in all images. North is up. Left panel (M = 0): the "signal" of structural changes 
. Middle panel (M = 0.009): best overall flattening of the rings using the magniflcation 
method. Right panel (M = 0.025 and the image scale has been doubled): best overall 
flattening of the lobes using the magniflcation method. Note that the expanding rings are 
absent (prominent) in projection on the north (south) lobe implying that the lobes are largely 
opaque at 0.6//m. 
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Fig. 3. — Ratio images of CRL 2688 shown in polar coordinates. Left: unsharp-masked ratio 
image in Cartesian (upper) and polar coordinates (lower), and assuming M — 0. Top right: 
as lower left except for a radial translation of the 2002.94 image by O'.'O? prior to forming the 
ratio image. Lower right: as lower left after a magnification of the 2002.94 image by a factor 
of 1.009. Although the rings are well flattened throughout most of the lobes in the upper 
and lower images, faint "ghosts" of the rings are more prominent within the white box. 



Fig. 4. — Left: sketch of the 2D model for the circumstellar dusty environment of CRL2688. 
All components except the "polar caps" are constrained in their shape by imaging data. 
Each of the model components is described in Appendix A. Components in the sketch are 
not shown exactly to scale; for the list of their parameters see Table 4. Right: color rendition 
of the model images with a non-linear contrast, blue = 0.6yum, green = l.l/xm, red = 1.6/im. 
The emerging light from the nebular core at 4.5/im is overplotted in red. 
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Wavelength X (/U.m) 



Fig. 5. — A fit of tlie predicted SED witli observed data points before (open circles) and 
after (solid circles) corrections for foreground extinction of Ay = 0.88 (U06). 
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Wavelength (jm\) 

Fig. 6. — Cross section of dust grains used in our model. Solid lines are absorption, dash-dot 
lines are scattering cross sections. Grain radii are indicated in the figure. A distribution of 
grain sizes for the smallest grains (lower curves) is a power law of exponent -3.5. 
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A. Appendix 

We describe the details of the major structural components used in or derived from the 
model. The absorption and scattering properties of the dust shown in Fig. [6] were adopted 
for the model, as discussed for each component of the model below. 

Star: The central star is not visible in images due to large dust opacity (~ 55) along 
the line of sight. Nonetheless, its basic spectral type can be derived from the molecular lines 
in the spectrum measured in the scattered light from bipolar lobes. Using this approach. 



Crampton et al. (1975) classified the stellar spectra as of type F5 la, while Cohen & Kuhi 



(1977) argue that it is F2 I. Thus we initially fixed the stellar temperature as a 7000K 



black-body spectrum which is similar to the temperatures suggested by these spectral types. 

In the end we found that several combinations of stellar spectral types and adopted 
reddening properties of the polar caps provide good fits to the nebular colors and the 
SED provided that we fix the luminosity at 5500 L0. For example, a very satisfactory 
fit is obtained if we fix T^, at 7000K and impose internal reddening consistent with ~ 2 
magnitudes of extinction in the dense polar caps. However, we were also able to create a 
comparably successful model with a 4000K black-body star of about the same luminosity 
and no internal reddening. If we adopt the early-F-type spectral type of about L* ~ 5500 
Lq and the necessary internal reddening then the spectral type of the star is between Fl-Ib 
and F3-Ib. (If the distance is doubled then the spectral type becomes Fl-Iab and F3-Iab.) 

Central dust free cavity: F-type stars are far too hot to form dust in their atmospheres. 
So we assumed that the innermost dust is located at the perimeter of a (presumed) spherical 
cavity of radius Rin in units of stellar radii. The value of Rin is determined in part by the 
star's luminosity L,,,, where L,,, = 5500 ±20% is measured by fitting the SED of the composite 
nebula at the adopted distance of 420 pc. R^ is basically set by the requirement that the 
temperature along the walls of the cavity, Tmax, is less than condensation temperature 
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of the dust Tcond, although Tcond is not known unless the chemical nature of the dust is 
somehow determined. Unfortunately, the high disk opacity precludes near-IR measurements 
of the color temperature of the cavity which might at least constrain the adopted value 
of Tmax- For the same reasons the size and geometry of the cavity cannot be measured. 
Moreover, we found that suitable models could be developed over a substantial range of 
adopted values of Tmax- Ultimately we simply adopted Tmax = 1000 K in the midplane. It 
then follows that Rin = 46.4/2,,= 0.56 msec with a formal error of 10% based on luminosity 
uncertainties. (Given the many assumptions, the actual uncertainty is obviously larger.) 

Disk: Imaging reveals a dense, flattened, and modestly articulated dusty structure in 
the midplane of CRL2688. Its irregular shape is exposed by multiply scattered starlight 
from the lobes (see Fig. [T]) and supported by the morphologies of H2 and radio molecular 
images noted earlier. For modeling purposes we call this structure an oblate "disk". The 
SED slope between 0.35 and 3.6 cm points toward emission from large dust grains, with 



sizes of 0.5 cm or more (Jura et al. (2000)). On the other hand, the 3.8 and 4.7/im images 



by G02 showed that the dust grains are not larger than ~ 0.2/im. 

This mismatch is a result of inhomogeneous spatial distribution of grain sizes, where 
large grains probably exist closer to the star in the irregular core of the dusty structure, 
while smaller grains dominate the opacity in the outer edges of the disk. We simplify the 
disk structure using single-grain dust of size 0.4yum with a radial density distribution that 
varies as r", where r is the distance form the star and a is a free parameter. We adopted 
an oblate ellipsoidal disk shape with bipolar lobes carved out at polar angles < 15°. For 
simplicity we assume that the edges of the bipolar cones are abrupt, so that the shadow 
that they cast in the envelope is sharply defined. The semi-major axis of the disk is ~ 7" 
and the semi-minor axis is ~ 1'.'8. The disk's midplane optical depth is adequate to hide the 
star at A < 4.5/im. 
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Polar caps: Polar caps are a vital construct of the model that collimates the light that 
both illuminates the lobes and escapes to form the searchlights. From the Earth the line of 
sight to the caps traverses the dense disk, so the caps are not fully visible directly in HST 
images in F606W and F814W. (Parts of the caps may be seen through patchy extinction in 
the FllOW and F160W images.) 

We noticed that without the caps, the scattering model for illumination from a central 
star of T^, ^ 7000 K produces lobes that are too bright and too blue to fit the data. Thus 
the starlight is modulated by an unobserved and compact region of dense dust. The only 
way to fit the shape of spectral continuum in the visible range with a 7000K star is to redden 
and partially absorb the light that enters the lobes to match the illumination expected from 
an illumination source >4000K . In essence, the caps are dust density enhancements at the 
base of lobes and, so, may be parts of a single entity. However, the requisite dust grains 
differ in size from the larger ones, 0.4/im used in the lobes and the disk. The dust in the 
caps consists of particles from from 0.01 to O.OSyum that follow a radius"^'^ size distribution 
as shown in the lower curve of Fig. |6| Without observations the density of these grains is 
difficult to characterize. We adopted a constant density pcap within a cone segment where 
the cone height Reap and the opening angle 6'cap are free parameters. This high-density cone 
is added to the densities of the lobe and disk present in this region. 

Lobes: The bipolar lobes are the brightest of all features detected in our images of 
CRL2688. The (identical) lobes have sharp edges in HST images, which enables precise 
tracing of their geometry for distances larger than about I'.'S from the position of the central 
star. Our model is symmetric with respect to the midplane. We analytically describe the 
shape of each lobe as an axially truncated cone with apex at the nucleus plus half of an 
ellipsoid of semi-major axis = 7" attached to the base of the cone (Fig. |4]) They join where 
their mutual widths are 1'.'26. The cones have a total opening angle of 30°. The internal 
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densities follow a power law, where b is varied to reproduce the observed radial trends 
in lobe surface brightness. We adopted an isotropic dust scattering function. Dust grain 
size in lobes must be dominated by grains larger than about 0.1/xm to assure that the lobes 
remain bright as seen in near infrared images. For simplicity we adopted the same dust 
properties as in the disk. 

Although the lobes are directly illuminated by the central star, their surface 
brightnesses are almost fiat owing to multiple interior scatterings. However, the northern 
lobe is observed to be slightly brighter then the southern lobe. The only way how we 
can achieve this asymmetry is through an axial inclination of angle i and the assumption 
of internal reddening within the envelope. A value of i ~ 15° provides a good match to 
the lobe ratio for the properties of the envelope that we assumed. Smaller values of i 
provide slightly better fits to the projected appearance of the arcs illuminated within the 
searchlights. 

Based on their HST imaging, S98b proposed a model for searchlight beams where the 
starlight initially escapes through a clear hole on the polar axis near the nucleus. A shadow 
zone within the searchlight is the result of shadows cast by a pair of extended knots lying at 
the outer tips of the the lobes. The existence of the knots is motivated by features seen at 
the tips of the lobes seen in H2 and CO images. In principle, the shadow-casting knots can 
be located at any radius along the polar axis. We created an annular aperture by locating 
a cone of enhanced density within each lobe whose opening angle, 10° in total, to fit the 
shadow edges within the searchlights. The dust within this inner cone is like that in the 
disk and outer lobes. Its density varies as r"^ where c is determined by fitting the darkness 
of the shadows. 

Envelope: CRL 2688 shows a faint spherical nebulosity out to angular radii r ~ 40" 
from the star. We call this region the nebular "envelope" within which thin spherical shells 
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mimic the conspicuous co-centric rings or arcs of varying size and brightness. For simphcity 
and speed we fixed the envelope radius to be twice the size of the lobes, ~ 15". This is 
big enough to analyze the properties of the rings and to collect all the fiux at mid-IR 
wavelengths where the SED reaches its maximum and begins to decline. 

We use dust of uniform radius in the envelope. The grain radius is a free parameter. 
The envelope density profile is constrained by its measured brightness in scattered light. 
To match the general properties of the observed envelope we assume that the dust density 
distribution oc (i.e., a mass loss rate that increases linearly with time for constant 
outfiow velocity, a common feature among pPNe Hrivnak & Bieging 2005| ). This produces 



a brightness profile ~ r~^, similar to the measured trend of ~ t--^.? fQ^j^^j Iqj gggb. The 
composite image on the right side of Fig. |4] shows a bluish color of envelope compared to 
the disk color, indicating that the envelope contains smaller grains than the disk and the 
lobes. However, the color of the searchlights is blue only in proximity to the blue lobes that 
illuminate them locally after multiple scattering. The color reddens at larger radii where 
starlight emerging from the caps is the primary light illumination source. This trend is not 
observed in our WFC3 images, perhaps because the grains in the envelope contain many 
smaller particles than those we used in the model. 

Shells: S98b concluded that the observed rings (represented by spherical shells in 
our model) are local density enhancements. They argued that a density enhancement by 
a factor of three yields a brightness increase of a factor of two. We used these density 
enhancement values relative to the ambient density of the background envelope that declines 
as r~^. Since observed arcs are broken and deformed, our aim is not to reproduce them in 
detail and we focus only to their general brightness properties. 

Spectral energy distribution: matching the SED is one of the greatest challenges — and 
successes — of our model. We fitted the SED after assuming that interstellar dust along the 
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line of sight reddens the visible and near-IR light, as described in §4. However, we note that 
we also managed to derive successful models of the SED without interstellar reddening, but 



this is merely of academic interest. We also plot centimeter fluxes from Jura et al. (2000) 
even though our model does not address the emission at radio wavelengths. Presumably the 
radio Amm and sub-mm flux is the direct result of thermal dust emission which is outside 
the scope of our model. 

Technical Details of the Model. The fitted model parameters are listed in Table 4. Some 
do not have units. The reason is that the model works with dimensionless variables, where 
the scale is set by the optical depth. Thus we do not parameterize the mass and density of 
dust particles in physical units. The absolute scale is set by the total radial optical depth 



in the midplane. A detailed discussion on this modeling issue is given in (Vinkovic 2003). 



The dust densities are not explicitly specified in the model because the absolute density 
values are scaled out in radiative transfer equations and replaced with scattering and 
absorption optical depths. The calculation determines the optical depth along an arbitrary 
radial line and an arbitrary wavelength. We use the optical depth ro.55 in the midplane 
at 0.55yum. The overall algorithm describing dust density at a point {x,y,z), for a given 
semi-major axis of disk, Rdisk (semi-minor axis = 0.25 Rdisk) and lobe Riobe (semi-minor 
axis = 0.18 Riobe), can be symbolically described as this: 

r=x^-|-y^+z^ 

if (r < 1) then density=0 

else if (x < z ■ tan(^Me/2) AND x^ < 0.18^ ■ {Rl^^^ - z^)) then { 

if (x < z ■ tan^O shadow f^) then density oc r'^ 

else density oc r'^ 
} else if (z^ < 0.25^- {Rjisj, - x^)) then density oc 
else { 
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density oar ^ 

if (r within shell) then density = 3- density 

} 

if (r < Reap AND x < z • tan(^cap/2)) density = density + Qcap 



Table 1: Log of HST/WFC3 Images of CRL2688 from G011580 



WFC3 




Filter center/ 


Pixel 


Exposure 


Archive 


Caiiiera 


Filter 


widtli (//111) 


Size n 


Time (s) 


Data Set 


UVIS1/SUB512 


F606W 


0.59/0.22 


0.040 


100 


IBlMOlOll 


UVIS1/SUB512 


F606W 


0.59/0.22 


0.040 


600 


IB1M01021 


UVIS1/SUB512 


F814W 


0.80/0.15 


0.040 


150 


IB1M01031 


UVIS1/SUB512 


F814W 


0.80/0.15 


0.040 


1000 


IB1M01041 


IR/SUB512 


FllOW 


1.15/0.44 


0.128 


70 


1B1M05RMQ 


IR/SUB512 


FllOW 


1.15/0.44 


0.128 


1178 


IB1M05020 


IR/SUB512 


F160W 


1.54/0.27 


0.128 


47 


IB1M05RNQ 


IR/SUB512 


F160W 


1.54/0.27 


0.128 


1178 


IB1M05010 
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Table 2: Log of HST F606W Images of CRL 2688 





Pixel 


Exposure 


Observation 


GO 


Archive 


Camera 


Size (") 


Time (s) 


Date 


Program 


Data Set 


WFPC2/WF 


0.10 


260 


07/17/95 


6221 


U2RC0303T 


WFPC2/WF 


0.10 


60 


04/09/98 


7423 


U4SHA101R 


WFPC2/WF 


0.10 


200 


04/09/98 


7423 


U4SHA102R 


WFPC2/WF 


0.10 


200 


04/09/98 


7423 


U4SHA103R 


WFPC2/WF 


0.10 


200 


04/09/98 


7423 


U4SHA104R 


WFPC2/WF 


0.10 


200 


04/09/98 


7423 


U4SHA105R 


ACS/WFC 


0.050 


110 


10/16/02 


9586 


J8GH55PLQ 


ACS/WFC 


0.050 


360 


10/16/02 


9586 


J8GH55PNQ 


ACS/HRC 


0.025 


200 


12/10/02 


9463 


J8DI81011 


ACS/HRC 


0.025 


200 


12/10/02 


9463 


J8DI81041 


ACS/HRC 


0.025 


200 


12/10/02 


9463 


J8DI81051 


WFC3/UVIS 


0.040 


100 


08/05/09 


11580 


IBlMOlOll 


WFC3/UVIS 


0.040 


600 


08/05/09 


11580 


IB1M01021 
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Table 3: Methods fc Results for Arcs &: Rings 



Cameras 


Interval 


Method 


M 


Exp. Age (y) 


WFPC2-WF3/WFC3 


2009.59 - 1995.54 


Difference 


0.013±0.003 


1076±30% 


WFPC2-WF3/WFC3 


2009.59 - 1998.27 


Difference 


0.013±0.003 


871±35% 


WFPC2-WF3/WFC3 


2009.59 - 1998.27 


Ratio 


0.012±0.003 


943±35% 


ACS-HRC/WFC3 


2009.59 - 2002.94 


Difference 


0.010±0.002 


665±20% 


ACS-HRC/WFC3 


2009.59 - 2002.94 


Ratio 


0.009±0.002 


739±10% 
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Table 4. List of model parameters 



Parameter description 


Parameter 


Value 


Comment 


lobe optical depth 
along polar axis 


'^lobes 


2 


required to explain the lobe colors 


disk optical depth 
in the midplane 


'^disk 


55 


it has to be large to hide the star in the near IR; additional 
constraint comes from the strength and spectral shape of 
mid IR fiux; variations within 10% are acceptable 


outer grid radius 


Rout 




arbitrarily fixed; the detected nebulosity extends further out 


stellar temperature 


T. 


7000K 


black body spectrum 


maximum dust 
temperature in 
the midplane 


T'max 


lOOOK 


poorly constrained; models with larger or smaller Tmax are 
allowed; this is maximum temperature of dust grains in 
the midplane 


disk density profile 




a = -0.5 


poorly constrained; variations ±100% are allowed 


disk semi-major axis 


-^disk 


300Ri„ 


this value of Rdisk corresponds to 7-8" 

(based on the luminosity derived from the BED fit) 


disk, lobe and shadow 
dust grain radius 




OAfim 


see Figure joj poorly constrained; 


lobe density profile 


Br'' 


B = 2x IQ-'^ 
b = 1 


constrained by lobes brightness; B has to be small enough 
to make lobes optically thin 


lobe semi-major axis 


Rlobe 


300i?i„ 


based on luminosity derived from the BED fit, this value 
of Rlobe. corresponds to 7-8" 


lobe opening angle 


Slobe 


30° 


constrained by lobe images 


shadow opening angle 


^shadow 


10° 


constrained by "searchlights" images 


shadow density 
profile 


Cr" 


C = 2000 
c = -3 


poorly constrained, but its optical depth has to be larger 
than 1 ; steep density profile puts the most of optical 
depth in the base of lobes 


polar cap density 


Qcap 


6500 


equivalent to the optical depth of ~2 at 0.55/xm, which 
dims and reddens the starlight entering the lobes 


polar cap radius 


Reap 




poorly constrained; optical depth of the cap is the key 
parameter 


polar cap opening 
angle 


^cap 


40° 


poorly constrained, but it must be > di^be lobe opening 
angle; optical depth of the cap is the key parameter 


envelope density 




D = 4.2 X 10" 


constrained by the envelope brightness profile 



profile 
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Table 4- 


-Continued 


Parameter description 


Parameter 


Value 


Comment 




envelope dust 
grain radius 




0.2fim 


see Figure joj poorly constrained; must be large enough 
to keep envelope brightness high in the near IR, 
but smaller than grains in the disk 


geometrical thickness 
of shells 




15Ri„(0'.'39) 


poorly constrained because arcs are of irregular shape 
and thickness 


radial distance 
between shells 


drshell 


35iJ,„(0"9) 


poorly constrained because arcs arc of irregular shape 
and thickness 


Derived parameters 


inclination angle 


i 


15° 


constrained by the brightness asymmetry between 
the northern and southern lobe 


bolometric flux 


Fbol 


1.0 ±0.2 X 10-8 


fit to the data gives 12% smaller flux, while this value 
takes into consideration the asymmetry in diffuse radiation 


luminosity (Lq) 




5500 ± 1100 


based on 420pc distance and SED fit 


stellar radius (-Rq) 


R* 


51 ± 5R, 
0.56 ± 0.05 msec 


based on luminosity, stellar temperature 
and distance = 420pc 


dust condensation 
radius 




46.4R*(0"026) 


critically depends on Tmax , after that on T0.55 
and optical properties of dust in disk 
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